Astronomy & Astrophysics manuscript no. Fontani'deuterationsurvey 


©ESO 2013 


January 18,2013 





Letter to the Editor 



Deuteration as an evolutionary tracer in massive-star formation * 

F. Fontani 1 ' 2 , Aina Palau 3 , P. Caselli 4 , A. Sanchez-Monge 5 , M. J. Butler 6 , J. C. Tan 6 ' 7 , 1. Jimenez-Serra 8 , G. Busquet 5 ' 

S. Leurini 10 , and M. Audard 1112 

1 ESO, Karl Schwarzschild Str. 2, 85748 Garching bei Munchen, Germany 

2 Institut de Radio-Astronomie Millimetrique, 300 rue de la Piscine, 38406 Saint Martin d'Heres, France 
e-mail: fontani@iram. fr 

3 Institut de Ciencies de l'Espai (CSIC-IEEC), Campus UAB-Facultat de Ciencies, Torre C5-parell 2, E-08193 Bellaterra, Spain 

4 School of Physics and Astronomy, E.C Stoner Building, University of Leeds, Leeds LS2 9JT 

5 Departament d'Astronomia i Meteorologia (IEEC-UB), Institut de Ciencies del Cosmos, Universitat de Barcelona, Marti i 
' Franques, 1, E-08028 Barcelona, Spain 

04 | 6 Department of Astronomy, University of Florida, Gainesville, FL 32611, USA 

5^ . 7 Department of Physics, University of Florida, Gainesville, FL 32611, USA 

■ 8 Harvard-Smithonian Center for Astrophysics, 60 Garden St., Cambridge MA 02138, USA 

9 INAF-Istituto di Fisica dello Spazio Interplanetario, Via Fosso del Cavaliere 100, 00133, Roma, Italy 

10 Max-Planck-Institut fur Radioastronomie, Auf dem Hugel 69, 53121, Bonn, Germany 
, 11 ISDC Data Center for Astrophysics, University of Geneva, Ch. d'Ecogia 16, 1290 Versoix, Switzerland 

■ 12 Geneva Observatory, University of Geneva, ch. des Maillettes 51, 1290 Versoix, Switzerland 



in 



Received date; accepted date 

ABSTRACT 



Context. Theory predicts, and observations confirm, that the column density ratio of a molecule containing D to its counterpart 
O || containing H can be used as an evolutionary tracer in the low-mass star formation process. 

Aims. Since it remains unclear if the high-mass star formation process is a scaled-up version of the low-mass one, we investigated 
whether the relation between deuteration and evolution can be applied to the high-mass regime. 

Methods. With the IRAM-30m telescope, we observed rotational transitions of N 2 D + and N 2 H + and derived the deuterated fraction 



o 



in 27 cores within massive star-forming regions understood to represent different evolutionary stages of the massive-star formation 
process. 

Results. Our results clearly indicate that the abundance of N 2 D + is higher at the pre-stellar/cluster stage, then drops during the 
formation of the protostellar object(s) as in the low-mass regime, remaining relatively constant during the ultra-compact HII region 
phase. The objects with the highest fractional abundance of N 2 D + are starless cores with properties very similar to typical pre-stellar 
\^ , cores of lower mass. The abundance of N 2 D + is lower in objects with higher gas temperatures as in the low-mass case but does not 

£f>l • seem to depend on gas turbulence. 

Conclusions. Our results indicate that the N 2 D + -to-N 2 H + column density ratio can be used as an evolutionary indicator both in 
ly-^ | low- and high-mass star formation, and that the physical conditions that influence the abundance of deuterated species likely evolve 

• , similarly along the processes that lead to the formation of both low- and high-mass stars. 

m . 

, Key words. Stars: formation - ISM: clouds - ISM: molecules - Radio lines: ISM 

J> ' 1 ■ Introduction Df rac and the core evolution in the low-mass case. This relation 
• *h ' predicts that Df rac increases when the starless core evolves to- 
^ I The study of deuterated molecules is an extremely useful probe wards the onset of grav it a tional collapse because, as the core 
»H ■ of the physical conditions in star-forming regions. Deuterated density profile becomes more and more centrally peaked, freeze- 
species are readily produced in molecular environments char- out of co increa ses in the core centre and hence the abun- 
acterised by low temperatures (T < 20 K) and CO depletion dance of deuterated molecules is greatly enhanced. When the 
(Millar et al. 1989). These physical/chemical properties are com- young stellar ob j ect f orme d at the core centre begins to heat its 
monly observed in low-mass pre-stellar cores (starless cores on sulTOUnd ings, the CO evaporated from dust grains starts to de- 
the verge of forming stars), where the deuterated fraction (here- stroy me deuterated species and D frac decreases. Observations 
after D frac ) of non-depleted molecules, defined as the column of both sta rless cores and cores with already formed protostars 
density ratio of one species containing deuterium to its counter- confirm the theoretical predictions in the low-mass regime: the 
part containing hydrogen, is orders of magnitude larger than the p re - S tellar cores closest to gravitational collapse have the high- 
[D/H] interstellar abundance (of the order of ~ 10 5 , Oliveira et est Dfrac ( Cra psi et al. 2005), while D frac is lower in cores associ- 
al. 2003). Caselli (2002a) found a theoretical relation between ated with Class 0/1 protostars, and the coldest (i.e. the youngest) 
objects possess the largest Z)f lac , again in agreement with the pre- 



* Based on observations carried out with the IRAM 30m telescope, dictions of chemical models (Emprechtinger et al. 2009, Friesen 
IRAM is supported by INSU/CNRS (France), MPG (Germany), and et al. l2010h . On the basis of these results, Df rac can be considered 
IGN (Spain). 
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as an evolutionary tracer of the low-mass star formation process 
before and after the formation of the protostellar object. 

Can this result be applied to the high-mass regime? This 
question is difficult to answer because the massive-star forma- 
tion process is still not well-understood: large distances (> 1 
kpc), high extinction and clustered environme nts mak e observa- 
tions of the process challenging (Beuther et al l2007al Zinnecker 
& Yorke 2007). Observationally, the study of Pillai et al. (2007), 
performed with the Effelsberg and IRAM-30m telescopes, mea- 
sured high values of Df rac (~ 0.2) from deuterated ammonia 
in infrared dark clouds, which are understood to represent the 
earliest stages of massive star and stellar cluster formation. In 
more evolved objects, from IRAM-30m observations, Fontani et 
al. (2006) measured smaller values of Df mc (~ 10~ 2 ) from the 
ratio N2D + /N2H + , which are nevertheless much larger than the 
D/H interstellar abundance. Despite these efforts, no systematic 
study of the [D/H] ratio across all stages of high-mass star for- 
mation has yet been carried out. 

In this letter, we present the first study of the relation be- 
tween deuterated fraction and evolution in a statistically signif- 
icant sample of cores embedded in high-mass star forming re- 
gions spanning a wide range of evolutionary stages, from high- 
mass starless core candidates (HMSCs) to high-mass protostellar 
objects (HMPOs) and ultracompact (UC) HII regions (for a def- 
inition of these stages see e.g. Beuther et al l2007al) . This goal 
was achieved by observing rotational transitions of N2H + and 
N2D + with the IRAM-30m telescope. We chose these species 
because N2D + can be formed from N2H + only in the gas phase, 
tracing cold and dense regions more precisely than deuterated 
N H3, wh ich can also be formed on dust grains (e.g. Aikawa et 
al. l2005h and then evaporates by heating from nearby active star- 
formation. Even though the observations are obtained with low 
angular resolution, the objects observed in the survey were care- 
fully selected to limit as much as possible any emission arising 
from adjacent objects. 

2. Source selection and observations 

The source list is in Table IA-11 where we give the source coor- 
dinates, the distance, the bolometric luminosity, and the refer- 
ence papers. We observed 27 molecular cores divided into: ten 
HMSCs, ten HMPOs, and seven UC HII regions. The source 
coordinates were centred towards either (interferometric ) in- 
frared/millimeter/centimeter continuum peaks or high-density 
gas tracer peaks (NH 3 with VLA, N 2 H + with CARMA or PdBI) 
identified in images with angular resolutions comparable to or 
better than 6", either from the literature or from observations 
not yet published. In general, we rejected objects whose emis- 
sion peaks were separated by less than ~ 8" from another peak 
of a dense molecular gas tracer. This selection criterion was 
adopted to avoid or limit as much as possible the presence of 
multiple cores within the IRAM-30m beam(s). The evolution- 
ary stage of each source was established based on a collection 
of evidence: HMSCs are massive cores embedded in infrared 
dark-clouds or other massive star forming regions not associated 
with indicators of ongoing star formation (embedded infrared 
sources, outflows, masers); HMPOs are associated with interfer- 
ometric powerful outflows, and/or infrared sources, and/or faint 
(S v at 3.6 cm < lmJy) radio continuum emission likely tracing a 
radio-jet; and UC HIIs must be associated with a stronger radio- 
continuum (S v at 3.6 cm > 1 mJy) that probably traces gas pho- 
toionised by a young massive star. We did not include evolved 
HII regions that have already dissipated the associated molecu- 
lar core. We also limited the sample to sources at distances of 



less than ~ 5 kpc. We stress that the three categories must be 
regarded with caution because it can be difficult to determine 
the relative evolutionary stage. This caveat applies especially to 
HMPOs and UC HII regions, whose evolutionary distinction is 
not always a clear cut (see e.g. Beuther et al. l2007al) . Among 
the HMSCs, three sources (AFGL5142- EC, 05358-mm3, and 
I22134-G) have been defined as "warm" in Table lA^TI we ex- 
plain the peculiarity of these sources in Sect. [3] The observations 
of the 27 cores listed in Table IA-11 were carried out with the 
IRAM-30m telescope in two main observing runs (February 2 
to 4, June 19 to 21, 2010), and several additional hours allocated 
during three Herapool weeks (December 2009, January 2010, 
and November 2010). We observed the N 2 H + (3-2), N 2 H + (1- 
0), and N2D + (2-1) transitions. The main observational parame- 
ters of these lines are given in Table IA-21 The observations were 
made in wobbler-switching mode. Pointing was checked every 
hour. The data were calibrated with the chopper wheel technique 
(see Kutner & Ulich 119811) . with a calibration uncertainty of 
~ 20 - 30%. The spectra were obtained in antenna temperature 
units, r?, and then converted to main beam brightness tempera- 
ture, T MB , via the relation T* A = T MB r/ MB , where rj MB = B eS /F eS 
is 0.74 for N 2 D + (2-1), 0.53 for N 2 H + (3-2) and 0.88 for N 2 H + 
(1-0) lines, respectively. All observed transitions possess hyper- 
fine structure. To take this into account, we fitted the lines us- 
ing METHOD HFS of the CLASS program, which is part of the 
GILDAS software-^ developed at the IRAM and the Observatoire 
de Grenoble. This method assumes that all the hyperfine compo- 
nents have the same excitation temperature and width, and that 
their separation is fixed to the laboratory value. The method also 
provides an estimate of the optical depth of the line, based on 
the intensity ratio of the different hyperfine components. For the 
faintest N2D + lines, for which the hfs method gives poor results, 
the lines were fitted assuming a Gaussian shape. 

3. Results and discussion: is deuteration an 

evolutionary indicator of massive star formation? 

The spectra of N 2 D + (2-1) and N 2 H + (3-2) for all sources de- 
tected in N2D + are shown in Figures IB- II - IB -6 1 We detected 
N2H + (3-2) emission in all sources. We also found a remarkably 
high detection rate in the N 2 D + (2-1) line: 100% in HMSCs, 
64% in HMPOs, and 100% in UC HII regions. Such a high de- 
tection rate indicates that deuterated gas is present at every stage 
of the massive star and star cluster formation process, even in 
the surroundings of UC HII regions where the gas is expected 
to be hotter and more chemically evolved. Even though for 12 
sources we also observed the N2H + (1-0) transition, we al- 
ways computed the column density of N2H + and the deuterated 
fraction from the (3-2) line given its smaller telescope beam, 
to limit the contribution of nearby sources as much as possi- 
ble. An overall presentation of the data obtained, and a deeper 
analysis of all physical parameters, will be given in a forth- 
coming paper. We derived the N2H + and N2D + column den- 
sities, A^(N2H + ) and jV(N2D + ), from the line integrated inten- 
sity f ollowing the method described in the appendix of Caselli 
et al. (I2002bl) . Thanks to the selection criteria for our sources, 
for which interferometric maps of dense gas are available for 
most of the regions, a first estimate of the filling factor could be 
computed. However, because maps of the two transitions used 
to derive Df mc have not yet been performed (except for 122 1 34- 
VLA1), the source size was determined from interferometric 

1 The GILDAS software is available at 
|http://ww.irarnT fr/IRAMFR/GILDAS 
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measurements of NH3(2,2). This assumption seems reasonable 
because this line traces gas with physical conditions similar to 
those of N 2 H + (3-2) and N 2 D + (2-1). To take into account the 
possible effects of the evolutionary stage on the source size, we 
also computed an average diameter for each evolutionary group. 
This turns out to be: 6.5" for HMSCs, 4.1" for HMPOs, and 
5.5" for UC HIIs (Busquet 2010, Busquet et al. 201 1, Sanchez- 
Monge 2011, Palau et al. 2007, 2010). We stress that these an- 
gular diameters are consistent with the (few) N 2 H + and N 2 D + 
interferometric observations published to date (e.g. see the case 
of IRAS 05345+3 157, Fontani et al. 2008). The N 2 H + and N 2 D + 
column densities, their ratio (Z)f rac ), as well as the line parame- 
ters used in the derivation of the column densities, are listed in 
Table |A^31 

The method assumes a constant excitation temperature, 
T ex . For the N 2 H + lines, T ex was derived directly from the 
parameters given by the hyperfine fitting procedure corrected 
for the filling factor (see the CLASS user manual for details: 
|http://kam .fr/IRAMFR/GILDAS/d o^/html/class-html/class.htmi| . 
The procedure, however, cannot provide good estimates for 
optically thin transitions or transitions with opacity (r) not 
well-constrained (e.g. with relative uncertainty larger than 
30%). For these, we were obliged to assume a value for r ex (for 
details, see the notes of Table fA-3l >. For the N 2 D + (2-1) lines we 
were unable to derive r ex from the fitting procedure for almost 
all sources because t is either too small or too uncertain. In 3 
cases only was the optical depth of the N 2 D + (2-1) transition 
well-determined, and so is T ex : in two of these objects we 
found a close agreement between the estimates derived from 
the N 2 D + (2-1) and the N 2 H + (3-2) transitions. Therefore, the 
N 2 D + column density of each source was computed assuming 
the same T ex as for N 2 H + . Since N 2 D + (2-1) and N 2 H + (3-2) 
have similar critical densities and we measure similar T ex for 
both transitions, the two lines approximately trace similar 
material, so that computing Df rac using them is a reasonable 
approach. The N 2 H + column densities are on average of the 
order of 10 13-14 cm -2 , and the N 2 D + column densities are of 



order 10 



12-13 



cm 2 . Both values are consistent with similar 



observations towards massive star forming regions (e.g. Fontani 
et al. 120061) . The measured T ex corrected for filling factor are 
between ~ 7 and ~ 50 K and agree, on average, with the kinetic 
temperatures measured from ammonia, except for the colder 
HMSCs for which they are a factor of ~ 2 lower. 

The deuterated fraction for the three evolutionary groups 
is shown in Fig. [1] where we plot A^(N 2 D + ) against A^(N 2 H + ). 
There is a statistically significant separation between the HMSC 
group, which has the highest average Df lac (mean value ~ 0.26, 
cr = 0.22), and the HMPOs and UC HII groups, which have sim- 
ilar average deuterated fraction: mean Z)f rac = 0.037 (cr = 0.017) 
for HMPOs, and mean D frac = 0.044 (cr = 0.024) for UC HII 
regions. Both are about an order of magnitude smaller than that 
associated with HMSCs. A closer inspection of the data using 
the Kolmogorov-Smirnov statistical test shows that the separa- 
tion in Df mc between the HMSC group and that including both 
HMPOs and UC HII regions is indeed statistically significant: 
the test shows that the probability of the distributions being the 
same is very low (P ~ 0.004). This is strong evidence that the 
two groups differ statistically. Therefore, massive cores without 
stars have larger abundances of N 2 D + than cores with already 
formed massive (proto-)stars or proto-clusters. The abundance 
of N 2 D + , however, seems to remain constant, within the uncer- 
tainties, after the formation of the protostellar object until the 
UC HII region phase. That Df mc is of the order of ~ 0.2 - 0.3, 
on average, in HMSCs, and then drops by an order of magnitude 




10' 10" 

N[N 2 H + ] (x 10 12 cm" 2 ) 

Fig. 1. N 2 D + column density versus N 2 H + column density. Blue 
symbols correspond to HMSCs (triangles: "warm" cores, see 
text); green squares show HMPOs (open squares are upper lim- 
its); black asterisks correspond to UC HII regions. The two lines 
indicate the average values of Df mc for the HMSC group (i.e. 
0.26) and that of both the HMPO and UC HII groups (i.e. 0.04). 
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Fig. 2. Deuterated fraction, £> frac = JV(N 2 D + ),W(N 2 H + ), as a 
function of several parameters: kinetic temperature (a), N(N 2 H + ) 
(b), N 2 H + (3-2) line width (c) and N 2 D + (2-1) line width (d). 
The symbols have the same meaning as in Fig. Q] For some 
sources, the errorbars are not visible because they are smaller 
than the symbol size. In panel (a), only the sources with temper- 
ature derived from VLA ammonia observations are plotted. 



after the onset of star formation, indicates that the physical con- 
ditions acting on the abundance of deuterated species (i.e. den- 
sity and temperature) evolve similarly along both the low - and 
high-mass star formatio n proc esses (see e.g. Crapsi et al. 120051 
and Emprechtinger et al. 120091) . 

Another interesting aspect emerging from Fig. [T]is that the 
three HMSCs defined as "warm" in Table lA4l(AFGL5142-EC, 
05358-mm3, and I22134-G, marked as triangles in the figure) 
have Df lac almost an order of magnitude smaller than the others. 
These differ from the rest of the sub-sample of HMSCs because 
they have temperatures 7\ > 20 K (see Table IA-31 and panel 
(a) in Fig. |2). High angular resolution studies indicate th at they 
could be extern ally he ated (Zhang et al. l2002i Busquet [2010, 
Sanchez-Monge 1201 lh . so that they are likely to be perturbed 
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by nearby star formation and we expect their properties to be 
different from those of the other, more quiescent cores. An an- 
ticorrelation between Df lac and the distance to heating sources 
such as embedded protostars was fo und in the cluster-forming 
Ophiuchus-B clump by Friesen et al (2010). Our study tends to 
confirm the Friesen et al.'s finding, even though the poor statis- 
tics does not allow us to drive firm conclusions. We also point 
out that the four cores selected from the Butler & Tan (2009) 
work (G034-G2, G034-F1, G034-F2, G028-C1) have the high- 
est values of all measured Df lac and lie in infrared-dark regions, 
away from active star formation. These four cores are hence very 
similar to the prototype low-mass 'pre-stellar cores' (e.g.L1544, 
L694-2, see Crapsi et al. 2005) and we propose that these are 
good 'massive pre-stellar core' candidates. 

In Fig. |2l we plot Df mc as a function of several parame- 
ters: the kinetic temperature, the N2H + column density, and the 
line widths derived from both N2H + and N2D + . To search for 
possible (anti-)correlations between these parameters, we per- 
formed two statistical tests: the Kendall's t and the Spearman's 
p rank correlation tests Q. For 7\, the tests were applied to all 
sources in our survey with gas temperature derived from VLA 
interferometric ammonia observations (see Table IA-31 . As it can 
be inferred from panel (a) of Fig. [2] Df lac and 7\ are slightly 
anti-correlated (r = -0.38, p- -0.50), and Df mc is also anti- 
correlated with the N2H + column density (t = -0.43, p- -0.60, 
panel (b) in Fig. [2j. We also find a very faint anticorrelation be- 
tween D frac and the N 2 FT line width (t = -0.17, p = -0.23) and 
between D iac and the N 2 D + line width (r = -0.25, p = -0.35) 
(panels (c) and (d) in in Fig. [2] respectively). In particular, this 
latter is difficult to trust being affected by large uncertainties in 
the N2D + line widths. Emprechtinger et al. (2009) suggested that 
in low-mass star forming cores the deuteration is higher in colder 
and more quiescent cores, according to the predictions of theo- 
retical models. A similar trend was found also in a small s ample 
of seven massive star-forming clumps by Fontani et al. (2006) 
including both HMPOs and UC HII regions but not HMSCs. 
That the warmer sources have lower Df rac is not surprising and 
can be explained by the CO freeze-out and the chemical reac- 
tions leading to the enhancement of deuterium abundance be- 
ing strongly depressed when the temperature increases (Caselli 
et al. 120081) . The lack of correlation between deuterated frac- 
tion and line widths tells us that the deuterium fractionation pro- 
cess is independent of the gas turbulence. This result agrees with 
high-angular r esolut ion observations of cluster-forming regions 
(Fontani et al. 2009, Busquet et al. 1201 Oh . but given the large un- 
certainties (especially on the N2D + line widths), the conclusions 
must be interpreted with caution. We speculate that the anticorre- 
lation between Df rac and /V(N2H + ) could indicate that, assuming 
that Df rac decreases in the protostellar phase, the N2H + column 
density increases during the younger and most embedded period 
of the protostellar phase, as suggested by Busquet (2010) for a 
different sample of sources. 

In summary, our findings indicate that the physical condi- 
tions acting on the abundance of deuterated species (i.e. den- 
sity and temperature) evolve similarly during both the low- and 
high-mass star formation process. To confirm this, several ques- 
tions however need to be answered: in HMSCs, do the N2D + 
and N2H + emission peak at dust emission peak as in low-mass 
pre-stellar cores? What is the nature of the N2D + emission in 
evolved objects (HMPOs and UC HII regions)? Is the emission 
extended or fragmented into several condensations (as found in 
the few massive star forming regions observed with interferom- 



2 (http://www.statsoft.com/textbook/nonparametric-statistics/ ) 



eters)? To answer these questions, higher angular resolution ob- 
servations are necessary. On the theoretical side, we also need 
to investigate this proposed evolutionary sequence using astro- 
chemical models. 
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Appendix A: Tables 

Table I A- 1 1 contains the list of the observed sources selected as 
explained in Sect. 2 of the main body text, and give some infor- 
mation extracted from the literature about the star forming re- 
gions in which the sources lie. Table lA^2l presents the observed 
transitions and some main technical observational parameters. 
Table I A-3l shows the results of the fitting procedure to the N2D + 
(2-1) and N2H + (3-2) lines (see Sect. 2 of the main body text) 
of all sources, and the physical parameters derived from these 
results, namely the N2H + and N2D + column densities and their 
ratio, Df mc . Other parameters discussed in Sect. 3 of the main 
body text are also listed. 



6 



Fontani et al.: Deuteration in massive-star formation 



Table A-l. list of the observed sources. Col. 4 shows the velocity at which we centred the spectra, corresponding to the systemic 
velocity. Cols. 5 and 6 give the source distance and bolometric luminosity (respectively) of the associated star forming region. This 
latter is a very rough first approximation of the core luminosity because it is based on infrared measurements having poor angular 
resolution. We adopt as source names those used in the reference papers listed in Col. 7. 



source 


RA(J2000) 


Dec(J2000) 


Vlsr 


d 




Ref. 




h m s 


o r n 


km s _1 


kpc 






HMSC 


I00117-MM2 fl 


00:14:26.3 


+64:28:28 


-36.3 


1.8 


io iA 


(1) 


AFGL5142-EC 01 " 


05:30:48.7 


+33:47:53 


-3.9 


1.8 


10 3.6 


(2) 


05358-mm3 B "' 


05:39:12.5 


+35:45:55 


-17.6 


1.8 


10 3.8 


(3,11) 


G034-G2(MM2) fl 


18:56:50.0 


+01:23:08 


+43.6 


2.9 


10 16r 


(4) 


G034-F1(MM8)" 


18:53:19.1 


+01:26:53 


+57.7 


3.7 


1Q 1.9r 


(4) 


G034-F2(MM7)" 


18:53:16.5 


+01:26:10 


+57.7 


3.7 




(4) 


G028-C1(MM9)" 


18:42:46.9 


-04:04:08 


+78.3 


5.0 




(4) 


I20293-WC 


20:31:10.7 


+40:03:28 


+6.3 


2.0 


10 3.6 


(5,6) 


I22134-G* 1 " 


22:15:10.5 


+58:48:59 


-18.3 


2.6 


10 41 


(7) 


I22134-B* 


22:15:05.8 


+58:48:59 


-18.3 


2.6 


io 41 


(7) 


HMPO 


I00117-MM1 


00:14:26.1 


+64:28:44 


-36.3 


1.8 


10 31 


(1) 


I04579-VLA1 C 


05:01:39.9 


+47:07:21 


-17.0 


2.5 


10 3.6 


(8) 


AFGL5142-MM 


05:30:48.0 


+33:47:54 


-3.9 


1.8 


10 3.6 


(2) 


05358-mml 


05:39:13.1 


+35:45:51 


-17.6 


1.8 


1Q 3.8 


(3) 


18089-1732 i 


18:11:51.4 


-17:31:28 


+32.7 


3.6 


1Q 4.5 


(9) 


18517+0437* 


18:54:14.2 


+04:41:41 


+43.7 


2.9 


10 41 


(10) 


G75-core a 


20:21:44.0 


+37:26:38 


+0.2 


3.8 


1Q 4.8 


(11,12) 


I20293-MMl fl 


20:31:12.8 


+40:03:23 


+6.3 


2.0 


10 3.6 


(5) 


121307° 


21:32:30.6 


+51:02:16 


-46.7 


3.2 


10 3.6 


(13) 


123385° 


23:40:54.5 


+61:10:28 


-50.5 


4.9 


io 4 - 2 


(14) 


UC HII 


G5. 89-0.39° 


18:00:30.5 


-24:04:01 


+9.0 


1.28 


10 V1 


(15,16) 


I19035-VLA1 


19:06:01.5 


+06:46:35 


+32.4 


2.2 


io 3 - 9 


(11) 


19410+2336° 


19:43:11.4 


+23:44:06 


+22.4 


2.1 


10 40 


(17) 


ONl° 


20:10:09.1 


+31:31:36 


+ 12.0 


2.5 


io 4 - 3 


(18,19) 


122134- VLAl fl 


22:15:09.2 


+58:49:08 


-18.3 


2.6 


10 41 


(11) 


23033+5951° 


23:05:24.6 


+60:08:09 


-53.0 


3.5 


10 40 


(17) 


NGC7538-IRS9 fl 


23:14:01.8 


+61:27:20 


-57.0 


2.8 


1Q 4.6 


(8) 



Notes. (fl) Observed in N 2 H + (3-2) and N 2 D + (2-1); {b) Observed in N 2 H + (1-0), N 2 H + (3-2), and N 2 D + (2-1); w Observed in N 2 H + (1-0 ) 
and N 2 D + (2—1 ); (m) "warm" HMSCs; ( r) Lum inosity of the core and n ot of the whole associated star-forming region (Rathborne et al. 2010); 
(1) Palau et al. {20 10) <2) Busquet et al. fcOllh (3) Beuther et al. d2007bf) (4 ) Butler & Tan d2009h (5) Palau et al. d2007l) (6) Busq uet et al. |2010) 
(7) Bu squet {201 Oh m Sanchez -Mong e et al. d2008h (9) Beuther et al d2004) <10) S chnee & Carpenter J2009t) (u) Sanchez-Monge d201lh (12) Ando 
et al. <f20TTb (13) Fontani et al. $M)4di) (14) Fontani et al. d2004bl) (15) Hunter et al. (2008) (16) Motogi et al. d2011h (17) Beuther et al. d2002h (18) Su et 
al. d2009h (19) Nagayama et al. d201lh 

Table A-2. Observed transitions and technical parameters 



molecular 


frequency 


HPBW 


Av" 


Bandwidth" 


transition 


(GHz) 


(") 


(kms 4 ) 


(km s" 1 ) 


N 2 H + (1-0) 


93.17376* 


26 


0.126 


230 


N 2 H + (3-2) 


279.5 lW 


9 


0.042 


77 


N 2 D + (2-1) 


154.21718 rf 


16 


0.076 


139 



Notes. (fl) Resolution (Av) and bandwidth of the spectrometer used (VESPA). <6) frequency of the main hyperfine component (F[F = 23 — > 1 2 
, Pagani et al. I2007T) (c) frequency of the F t F = 4 5 — » 3 4 hyperfine comp onent, having a relative intensity of 17.46%. (Crapsi et al. 2005) 
<d> frequency of the main hyperfine component (F\F = 2 3 — » 1 2, Pagani et al. l2007h 
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Appendix B: Spectra 

In this appendix, all spectra of N2D + (2-1) and N2H + (3-2) tran- 
sitions, for the sources detected in both transitions, are shown. 
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HMPO*; 10(1117 MM 1 


S S(0 1 1 


SO 89(0 on 


1 S0(0 041 


2.4(0.3) 


19.3(0.4)" 


< 0.07 






1 26(0 1) 


< 0.5 


< 0.04 


20" 


104S79 VI A 1 f*^t 


] 'j'tm 091 


•19 M (0 Ofil 

— J Z.J 1 ^U.UU ) 


1 8(0 2) 


<0 1 


74 rf 


< 14 






4(0 0S1 


< 1 5 


< 4 


74' 


AFPil S149-MM 


27 2(0 1) 


1 7 470(0 0041 

— 1 / .*-r / \J\\J .yjyj'-T J 


9 870(0 0061 


3 4(0 1 1 


^9 94(0 0] \ h 


40(0 071 


-17 7(0 1) 


2.3(0.3) 


6 1 (0 6) 


3.0(0.5) 


049(0 0091 


34 /f 


05358-mml 


30.3(0.1) 


-30.588(0.003) 


2.040(0.003) 


sen 


46(?l ft 


0.15(0.02) 


-30.7(0.3) 


2.5(0.6) 


7.2(0.7) 


1.2(0.2) 


0.017(0.004) 


39 /! 


1 SORQ 1 7^9 


3 1 4(0 2) 


1 7 40(0 (\C\A\ 


4 960(0 001 1 


<0 1 




79(0 CWt 


1 X ^(0 61 


^(li 


7 0(0 7) 


2 2(0 3) 


rn i (O oofti 


18^ 


1 SSI 7+0437 


19.2(0.1) 


99 44C0 031 


3 08(0 011 


<o. l 


47'' 


1 S(0 071 
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^( 1 1 
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096(0 0061 

\/A/Z,U^V/A/V7U J 




H7S cnrp 
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1 A 90(0 091 


i. oqcn 041 


<0 1 


96'' 


< 05 






2 4(0 2) 


< 69 


< 029 


96'' 


190993-MM1 


9Q ficn 1 1 


8 31 CO 021 


1 610(0 0091 


6 4(0 1) 


50. 1(0. 


64(0 031 


9 3 1 (0 081 
— y .3 1 iu.uo ) 


1 8(0 3) 


7 3(0 71 


5.4(0.8) 


07(0 011 


37'' 


121307 


6 5(0 1) 


__f.\ 9Q/n 091 


9 Qfi(0 0S1 


<0 1 


21'' 


< 06 






1 S(0 11 


< 4 


< 03 


2 j ' 


123385 


9.41(0.08) 


-64.88(0.01) 


3.02(0.03) 


<(). 1 


37 ,/ 


0.08(0.05) 


-63(31 


1.2(0.9) 


2 1 (0 21 


0.6(0.1) 


0.028(0.009) 


37 l 


I 1C HI! OS 80 39 


84. 1(0.3) 


S 870(0 0031 


S 900(0 0091 


<0. 1 


20 1 


S0(0 031 


7 4(0 1 1 


1 S(0 31 


1 9(91 


3 4(0 SI 


01 8(0 0031 

U.Ul O^W. \J\.U J 




1 1 903S-VI A 1 


13 7(0 2) 


1 7 69(0 on 


4 49(0 031 


<0.1 


1 9 0(0 ?1# 


1 8(0 031 


17 5(0 1) 


1 4(0 31 


3.2(0.3) 


1 3(0 31 


04(0 01 1 


39* 


1 9410+7336 


33.3(0.1) 


7 90(0 00 SI 


3 SS0(0 001 1 

D —j J>\}\\.} .\j\3 y ) 


<(). 1 


21'' 


S9(0 0S1 


7 S(0 081 


1 S(0 31 


7 4(0 71 


3.5(0.5) 


047(0 0091 


21* 


ONI 


30 QCfl 71 


9 49 U0 0011 


A 340(0 001 1 


<0.1 


26 ,/ 


17(0 0S1 


3 3(0 91 


3. 13(1.65) 


6 6(0 7) 


1 1 (0 2) 


01 7(0 0041 


26 7 " 


19? 1 "KA- VI A 1 

1 Z.Z, 1 JT - V L-rt. 1 


9 S(0 081 

Z,.-J ^17 A70 ^ 


33 1 1 (0 011 


1 30(0 061 


1 9(0 41 


1 6(0 Sl^ 


08(0 Oil 


34.2(0.4) 


7(1 0) 


1 0(0 1) 


8(0 91 


08(0 071 

U.UO^U.UZ f 


47 A 


73033+S9S1 

ZJUJJTJ7J 1 


19.6(0.1) 


67 77(0 031 


3 S30(0 0021 


<0. 1 


25^ 


S9(0 071 


68 9(0 SI 


8(0 SI 


4 4(0 41 


3.5(0.5) 


08(0 071 

U.UO^U.UZ f 


25* 


NGC7538-IRS9 


13.9(0.1) 


-71.814(0.006) 


3.42(0.03) 


<0J 


20^ 


0.14(0.02) 


-72(2) 


2.6(0.5) 


3.2(0.3) 


0.9(0.2) 


0.030(0.008) 





Not es. (+) r e x associated with uncertainties are calculated from the data, the others are assumed as explained in the footnotes; (lt) 7\ are derived from ammonia rotation temperatures following Tafalla 
et al.|2004; (fl) r ex ~l/2 T k , based on the results derived in this work for the HMSCs with well constrained T ex and 7\, except the "warm" ones; ( ' derived from the (3-2) line parameters obtained from 
the hyperfine fitting procedure described in Sect. |2]and corrected for filling factor (see also the CLASS User Manual (http://iram.fr/IRAMFR/GILDAS/doc/html/class-html/class.htmll; (c) average 
value for HMSCs; (< " T ex ~ T k , based on the results derived in this work for the HMPOs with well constrained T cx and 7\, and on the findings of Fontani et al. 12006) towards a sample of m assive 
cores containing both HMPOs and UC HII regions; ^ average value for HMPOs; ( " average value for UCHIIs; {g) computed from the N 2 H + (1-0) line; from VLA observations (Busquet l2010l 
San chez-M onge l201 lL Palau personal comm.); (,) from Effelsberg observations (Pillai et al. 2006); {j) from Effelsberg observations (Molinari et al. 1996); <k) from Effelsberg observations (Sridharan 
et al l2002h ; (,) from VLA observations (Fontani et al. 2004b); < " l) from Effelsberg observations (Jijina et al. 1999); (T, r) N 2 H + parameters are derived from the (1-0) line, for which Jr MB di; is the 
integrated intensity of the isolated component (F t F = 1,0-1,1) appropriately normalised and assumed to be optically thin. However, because observations were obtained under very bad weather 
conditions, the £) fl!U . upper limit derived from this source is unreliable and we decided to exclude it in the following analysis. 
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HMSCs 

N 2 D (2-1) N 2 H + (3-2) 




V LSR (km s 1 ) 

Fig. B-l. Spectra of N2D + (2-1) and N2H + (3-2) obtained towards the sources classified as HMSCs. For each spectrum, the velocity 
interval shown is ±10 km s from the systemic velocity listed in Table lA^TI The y-axis is in main beam brightness temperature 
units. In the spectra of 1001 17-MM2, the vertical bars show the position of the hyperfine components. 
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HMSCs 

N 2 D (2-1) N 2 H + (3-2) 




V LSR (km s 1 ) 

Fig. B-2. Fig. EHJcontinued. 
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HMPOs 

N 2 D (2-1) N 2 H + (3-2) 




V LSR (km s 1 ) 

Fig.B-3. Same as Fig HO for HMPOs. 
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HMPOs 

N 2 D (2-1) N 




Fig. B-4. Fig. IB-3l continued. 
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UC Hlls 

N 2 D (2-1) N 2 H + (3-2) 




V LSR (km s 1 ) 

;. B-5. Same as Fig. IB-ll for UC HII regions. 
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Fig. B-6. Fig. IB-5l continued 
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